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Abstract. We study the electronic spectra of commensurate and incommen- 
surate double-wall carbon nanotubes (DWNTs) of finite length. The coupling 
between nanotube shells is taken into account as an intershell electron tunneling. 
Selection rules for the intershell coupling are derived. Due to the finite size of 
the system, these rules do not represent exact conservation of the crystal momen- 
tum, but only an approximate one; therefore the coupling between longitudinal 
momentum states in incommensurate DWNTs becomes possible. The use of the 
selection rules allows a fast and efficient calculation of the electronic spectrum. 
In the presence of a magnetic field parallel to the DWNT axis we find spectrum 
modulations which depend on the chiralities of the shells. 
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1. Introduction 

Due to their unusual physical properties, cf. e.g. [U [2], carbon nanotubes have 
become promising building blocks for nanotechnology applications and have attracted 
a lot of attention since their discovery. Carbon nanotubes can be single-walled 
(SWNT) or multi- walled (MWNT), depending on whether they consist of a single 
or of several graphcne sheets wrapped onto coaxial cylinders (so called "shells"), 
respectively. Electronic properties of SWNTs have been mostly understood pQ. 
For example, SWNTs are usually ballistic conductors [2], and whether a SWNT 
is metallic or semiconducting is solely determined by its geometry. However, the 
situation is much less clear for MWNTs. In fact, due to the additional shells, MWNTs 
exhibit qualitatively different properties than SWNTs. Except for few experiments, 
see e.g. @JH], MWNTs are typically diffusive conductors [Sill], with current being 
carried by the outermost shell at low bias [7J [8] and also by inner shells at high 
bias [9]. A recent experiment reported that the intershell conductance is quite weak 
and consistent with the tunneling through the orbitals of nearby shells [TU]. The 
difficulty in a theoretical description of MWNTs lies in the fact that the coaxial shells 
have usually different chiralities. In such case MWNTs are intrinsically aperiodic, 
since a common unit cell for the whole object cannot be defined due to the respective 
symmetries of individual shells. 

The simplest system in which the inter-shell effects can be studied is a double-walled 
nanotube (DWNT). The two shells are coupled by weak van der Waals interactions, 
which give rise to an inter-shell electron tunneling. DWNTs have been studied in 
various approaches. By using ab-initio methods on graphite, effective inter-layer 
hopping integrals have been found [11] , closely matching experimental results jT2l [13] . 
In the calculations involving nanotubes the hopping parameters are usually considered 
to be similar to those in graphite. An ab-initio study of multiwall nanotubes [14] 
correctly predicted the intershell distance and freedom of telescopic and rotational 
motion of the shells, later confirmed experimentally [15]. Commensurate DWNTs have 
been thoroughly analyzed, and their electronic spectra [HI [17j [HI [19] and transport 
properties [20] have been discussed. Some authors investigated also the properties 
of incommensurate DWNTs, looking at spectral correlations [5TJ [55] and transport 
properties [H [H [H [25l |2B] or simulating their STM images [27] . Each of the 
transport studies refers in some way to the selection rules for the inter-shell coupling. 
They are mentioned in passing in [22j when discussing the tunneling between states in 
the inner and outer shell at Fermi points K and K' . The general analysis of MWNT 
conduction presented in [24] relies on the conservation of quasi-crystal momentum to 
prove that the conductance of a long MWNT is dominated by the outermost shell. 
In [25] the authors consider a long DWNT and calculate the intershell resistance, as 
coming only from the Coulomb drag, i.e. neglecting the inter-shell tunneling. They 
find selection rules for the coupling between momentum states in different shells. 
Since the interest in the above works is focused on the conduction, they explore the 
consequences of those rules mainly close to the Fermi level. 

When a uniform magnetic field is applied to a system, interesting and subtle effects 
occur, depending on the geometry and topology of the system, due to the Peierls 
phase [28j acquired by the electronic wavefunction. For electrons moving in spatially 
periodic potentials, if the flux through the elementary cell contains an irrational 
number of flux quanta, the periodicity is destroyed and the spectrum becomes fractal 
[29] . When the field is applied parallel to the axis of symmetry of a ring or cylinder, 
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it causes the Aharonov-Bohm effect or persistent currents [3U1 IS]- In nanotubes, 
because of their unique dispersion relation, the field can induce e.g. a periodic 
metal-semiconductor transition, predicted in 32J and observed in many experiments 
[531 |3"?I [551 [5tj] , The effects of a uniform magnetic field on the spectrum of a 
commensurate DWNT in the vicinity of the Fermi level have also been studied, taking 
into account several rotational configurations of the two shells [37] [38] . The tunneling 
coupling between shells of a DWNT modifies the spectra of the individual shells, 
introducing numerous avoided crossings, which in turn result in the depletion of the 
density of states (DOS) in one or more regions of the spectrum [35] • I n small fields 
this region lies close to the bottom of the valence band, but when the magnetic field 
increases, the influence of the intershell coupling is visible in higher energy ranges. 
In this work we extend to finite size DWNTs an approach presented in [23] in which the 
Hamiltonian of DWNTs is analyzed in the reciprocal space. We find the selection rules 
for the coupling between momentum states and estimate the amplitude of the coupling. 
This method has the advantage of being computationally fast, due to the action of 
the selection rules and can be applied to commensurate as well as incommensurate 
DWNTs. For short DWNTs both our method and the direct diagonalization of the 
tight-binding Hamiltonian in the real space yield spectra with the same positions 
of the van Hove peaks with some mismatches in their heights. For commensurate 
DWNTs in parallel magnetic field our result matches the results of 39J, where a 
similar system has been studied in the real space. We also calculate the electronic 
spectra in changing magnetic field for incommensurate DWNTs. We find a periodic 
closing and opening of the gap at the Fermi level, as well as a region with depleted 
density of states (DOS). This region evolves with the magnetic field in a complex way, 
determined by the geometry of the two shells. 

This paper is organized as follows. In section[2]we introduce various quantities needed 
for the characterization of the real and reciprocal space of graphene and of nanotubes. 
The intershell tunneling in DWNTs is studied in section [3j where the reciprocal space 
formula for the tunneling coupling is derived and analyzed. The coupling changes 
when a uniform magnetic field is applied and its influence on the energy spectrum is 
studied in sectional Section [5] concludes this work. 

2. Direct and reciprocal lattice structure of DWNTs 

2.1. Graphene 

Various nanostructures, such as nanotubes, graphene ribbons or nanocones [40] [41] 
can be treated as fragments of a graphene sheet (figure [1]) with appropriate boundary 
conditions. For later purposes we briefly recapitulate how to characterize the graphene 
lattice and its electronic spectrum and how to adapt this description to the case of 
carbon nanotubes. The honeycomb lattice of graphene is generated by two basis 
vectors of equal length, and the angle between them is 60°. We choose their Cartesian 
coordinates as 



where ao = 1.42A is the length of a C — C bond. The elementary cell contains 
two atoms which generate the two sublattices of graphene through the translations 
by multiples of 0,1,0,2- The atoms A and B in the elementary cell are shifted with 




(1) 
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respect to the origin of coordinates by vectors T^.Tg, respectively. We will refer 
to these vectors as the sublattice shifts. In the Cartesian coordinates which we have 
chosen, they are given by 

T A = (0,0), TB = (0,Oo). (2) 

The generators of the reciprocal lattice are 

b - ( 2n 2lT \ b - ( 4?r ^ (3) 

V V3a 3a / \ ' 3a / 

A common starting point for the calculation of the band structure of graphene is the 
tight-binding model for nonintcracting p z electrons pQ , described by the Hamiltonian 

H = ^2 lQ c\ a c ja , (4) 

(if 

where 70 ~ — 2.7eV is the hopping integral in graphene; i,j are the p z orbitals of 
carbon atoms at positions i and j, respectively; a denotes the electron spin and the 
sum runs over nearest neighbours in the real space. The dispersion relation can be 
derived (see section [3] for details). It reads 



e v (k) = 2/7oy3 + 2 cos(fc • a{) + 2 cos(fc • a 2 ) + 2 cos(fe • (a 2 — ai)), (5) 

where v — +1 in the conduction band and v = — 1 in the valence band. This dispersion 
relation has the characteristic shape of a double crown, with six Fermi points - only 
two of them being geometrically inequivalent. A fragment of the atomic lattice of 
graphene and its reciprocal lattice are shown in figure [TJ 




Figure 1. Direct and reciprocal lattice of graphene. a) Atomic structure of the 
honeycomb lattice with two sublattices A and B and the lattice generators a\ , ai . 
Vectors d{ connect the atoms from sublattice A with their nearest neighbours, b) 
The first Brillouin zone of graphene and the reciprocal lattice generators 61 , 62 . 
The background is a greyscale map of the negative part of the dispersion relation 
in 70 units. 



2.2. Single-wall nanotube (SWNT) 

A single-wall nanotube can be described as a rectangular patch of graphene with two 
opposite sides joined together by periodic boundary conditions (figured). The vector 
defining the circumference of the SWNT is called chiral vector and is uniquely defined 
by two coordinates in the basis of lattice generators 

Ch = miOi + m 2 a 2 , mi,m 2 eZ. (6) 
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Because of the hexagonal symmetry of graphene this notation is redundant. In 
particular, the nanotube with (—mi, — 777,2) is identical to the one with (7711,7712), and 
(7772, 777!) is its mirror image. The convention is to keep mj > 7772 and 7772 > 0. In most 
nanotubes a chiral arrangement of atoms can be observed along the nanotube. There 
are only two combinations of parameters which describe achiral nanotubes: (m, 0) 
corresponding to so-called zigzag tubes and (777,777) corresponding to armchair tubes. 
The nanotubes can also be viewed as objects created by a repeated translation of a 
unit cell, defined by the vectors Ch and T (see figured): 

7771 + 2?772 , 27771 + 777 2 

/ = - aiH a 2 . (7) 

a-R a R 

Here oIr is the greatest common divisor of (mi + 27772) and (2mi + 7772). 
The boundary conditions around the circumference of the nanotube (in transverse 
direction) are always periodic {PBC), There are two ways of dealing with the boundary 
conditions along the nanotube axis (in the longitudinal direction), resulting in the 
same spectrum. One way is to consider open boundary conditions (OBC) with the 
wavefunctions defined on the length of the nanotube. The other way is to consider 
periodic boundary conditions on a nanotube twice that length and to choose only 
the energy eigenfunctions which are antisymmetric with respect to the center of the 
extended tube. Physically it means that we choose only those wavefunctions which 
are reflected from the end of the original tube (or the center of the extended tube) 
with opposite phase. This restriction removes both the level degeneracy caused by 
PBC and the fcii = eigenstate, which is symmetric with respect to the center of the 
extended tube. 

The boundary conditions cause the quantization of momentum 

k = (k x ,k ll ) = (^-l ± ,^, l±,l {l ez, (8) 

where L = M\T\ is the length of the nanotube, equal to M unit cells. Note that 
in ([5]) open boundary conditions along the nanotube axis have been assumed. In 
infinite nanotubes fc|| is continuous and the allowed momentum states are a set of 
lines of constant k± . Instead of working in the quantized hexagonal Brillouin zone 
of graphene, it is more comfortable to define a rectangular unit cell of the reciprocal 
space, with the area equal to that of the Brillouin zone and yielding the same energy 
spectrum (figure^)- We shall refer to it as the reciprocal cell. It is spanned by vectors 
b± and b|| given in the basis of graphene reciprocal lattice generators by 

27771 + 7772 , , 7771 + 2t77 2 , , n , 

b± = bi H b 2 , (9a) 

h =-^ 1 + ^6 2 . (96) 

Notice that the coordinates of b± are integer, therefore b± is always a reciprocal lattice 
vector. Note also that, since l± = k±R, the angular momentum is Kl±. 
The projection of the lines of constant k± on the dispersion relation reduces the full 
2D spectrum to a set of ID subbands, numbered by their value of angular momentum 
quantum number l± (figure^). The number S of subbands in one band, equal to the 
number of allowed values of is the number of graphene unit cells in the unit cell 
of the nanotube 

, , 2(mj + 77717772 + 777^) 
S(777l,777 2 ) = ~ . (10) 

a-R 
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Each subband has a positive and a negative energy branch, accounting for the presence 
of two atoms in the graphene unit cell. In the reciprocal cell all subbands contain equal 
number of fcy states. In finite nanotubes the ID subbands are further discretized, and 
a nanotube containing M unit cells has M longitudinal momentum values in each of 
the S subbands. Therefore < In < M and the allowed range of l± is [—5/2, 5/2). 




Figure 2. Characterization of an armchair nanotube. a) Unrolled nanotube 
patch (light grey) on a graphene lattice - the chiral vector is (5,5) and the 
nanotube has only four unit cells. The area of the unit cell, spanned by vectors 
Ch and T is marked in dark grey, b) The Brillouin zone of graphene (dashed 
lines) and the reciprocal cell (solid lines) of an infinite (5,5) tube with the allowed 
momentum states, c) The electronic subbands of an infinite (5,5) nanotube. 
Quantum numbers of the subbands in the conduction band E > are the same 
as their equivalents in the valence band. 

2.3. Double-wall nanotube (DWNT) 

A double- wall nanotube consists of two coaxial single- wall nanotubes, called also shells. 
The inter-shell distance A is typically of the order of 3.4A [42]. The coupling between 
two shells can be taken into account as an inter-shell tunneling of electrons. The 
implications of this tunneling will be explored in section [3] A schematic picture of a 
DWNT and its system of coordinates is shown in figure [3) 
In the present paper we will be using several systems of coordinates, each of them 




Figure 3. Systems of coordinates used to describe a location on a double-walled 
nanotube. a) Two graphene layers. The vectors p a , p b describing the relative 
position of the graphene patches of the two shells (armchair a and zigzag b) 
and shifts T a , T b between A and B sublattices in both shells are indicated, b) 
Schematic view of a DWNT and its system of coordinates (x±,Xt\). 
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suitable for a particular purpose. For the first three systems of coordinates we start 
from two graphene layers separated by a distance \Rf, — R a \. Each point on one of the 
constituent 2D graphene layers of a DWNT can be described either by the Cartesian 
coordinates (x, y) or by the nanotube patch coordinates (x±, sen). The third possibility 
is the system defined by (01,02), but this one is used only in the definition of the 
nanotube chirality Cn- When the nanotube is rolled, it becomes a 3D object and 
the most natural coordinate system is the cylindrical one. The cylindrical coordinates 
(r, ip, z) of a point on the shell {3 are related to the 2D nanotube coordinates (x±, Xu) 

by 

X I ft 

(r,<p,z)p = {Rp,~s —,x\\p). (11) 

In the reciprocal space we use only the 2D coordinates of the graphene layers. The 
vectors G in the reciprocal space can be expressed as (ni,na) in the graphene basis 
of (6i,6 2 ), i-e. G = n\b\ +n 2 6 2 , or as (G±,G\\) in the basis of (6j_,6||) spanning the 
reciprocal cell of a nanotube (see (|9a|) and (|9&|) ). 



3. Effective intershell coupling in DWNTs 

The starting point for our investigation of the consequences of the inter-shell electron 
tunneling is a tight-binding model for noninteracting p z electrons on each shell of the 
carbon nanotube pQ . The tight-binding Hamiltonian of a DWNT is obtained from that 
of two graphene sheets a and b placed on top of each other at a distance \R a — Rb\ by 
imposing periodic boundary conditions along the directions determined by the chiral 
vectors C/ m = (rni a ,m 2a ) and Chb = (mib>m2b) (see (j6])). The DWNT Hamiltonian 
is 

H = H + H t =J2J2 lv c li<j c pjo + ^«.^i c L-%- CT + H.c, (12) 

0a (ij) ija 

where the operators c^ jo . and c^ J(T are creation and annihilation operators of an 
electron with spin a on shell (3 at site j, respectively. Here (3 = a, b is the shell 
index and, as in (|4]), (ij) is a sum over nearest neighbors and 70 ~ — 2.7eV is the 
intrashcll nearest neighbor coupling. The spin-independent intershell coupling t i,6j is 
assumed to depend exponentially on the distance between two atoms, d(r a i,rf,j), as 

=tocos%e-^^)- A )/°*, (13) 

where to ~ — 0.34eV, A ~ 0.34 nm, 6ij is the angle between the p z -orbitals of the 
two atoms, and at ~ 0.45A [1] is a parameter controlling the range of the tunneling. 
We adopt here the second approach to the boundary conditions along the nanotube 
axis, described in section 12.21 We extend our DWNT to twice its original length, 
assume periodic boundary conditions, and reject all solutions which are symmetric 
with respect to the center of the extended nanotube. The sum over i,j runs therefore 
over the extended nanotube. 

It is convenient to express the Hamiltonian in the basis of plane waves in each 
individual shell [3 [33] . We introduce the electron operators 

*■ \ " ife-r-j t _ 1 \ " -ik-Tj t 
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where p = A, B is the index for the two interpenetrating sublattices in a graphenc 
sheet, and Np is the number of graphene unit cells on shell [3. The extended tube has 
twice as many atoms as the original one, hence the \[2 in the normalization factor. 
The Hamiltonian takes the form 



H 



^ 7fc C /3pfecr C /3p'fc(T 

fipka 



E 

k a k b p a p b <7 



1~PuPb(ka,kb)c 



t „ 
ap a k a a bp b k b a 



H.c. 



(14) 



where the mirashell coupling is 7^ = J2j=i lo e%k dj > with dj the vectors connecting 
an A sublattice atom to its three nearest neighbours in sublattice B (figure Hk) . The 
position of each atom in the graphene patch can be expressed as rp = R + Xp, with 
R a graphene lattice vector, Xp = pp + Tp p , where p a — p b describes the relative 
position of the two shells and Tp p is the appropriate sublattice shift, cf. figure [3] The 
elements of the intershell 2x2 coupling matrix can be expressed as [43] 



%>a.Pb(ka, kb) — 



AG a -X a -iG b -X 



G a G b 



Here G is the graphene reciprocal lattice vector G 
intershell coupling has the form 

1 



b t ka +G a M b +G b ■ (15) 

n\bi + n 2 b 2 = (n 1 ,n 2 ). The 



t 



la.Qb 



J dr a dr b e l{q >' rb - q « r " ) t ra 



Ai eU vmM' """ "' ""■ (16) 

with -Aceii the area of a graphene unit cell and the integral taken over the area of the 
system, in our case over the extended nanotube. For the purpose of calculating the 
energy spectrum, it is better to use the basis of the eigenstates (Bloch states) of the 
Hamiltonian pT]) in the absence of intershell coupling. This can be achieved by the 
unitary transformation 



7fc 

It* I 
1 



7fc 
7fcl 
1 



(17) 



The tunneling matrix elements between two Bloch states in different shells can be 
obtained as 

(18) 

Here v — =p is the index for two graphene bands corresponding to negative/positive 
energies ep jV (k) with [3 — a,b, where the dispersion relation of these bands is, cf. (J5j) , 



ep. v (k) = i>7o y3 + 2 cos(fc • ai) + 2 cos(A; • a 2 ) + 2 cos(fc • (02 
The electronic momenta are quantized according to the boundary conditions 



ax))- 



2tt 2ir 
kp ■ C h p = jq—lxp, kp ■ 2Lp = — — hp, l±p,hpeZ. (19) 

In order to calculate the inter-shell coupling (fT6]) we shall use nanotube coordinates, 
(R, xx/R, xn). The distance between two atoms a and b with cylindrical coordinates 
(R a , <p a ,z a ) and (R b , ip b , z b ) is thus 



d(r ai r b ) 



d(^ 
V Rb 



X_La 
Ra 



;X\\b 



\R a - Rb? +4i? a i? b sin 2 



X±b 

Rb 



X±g 
Ra 



+ (X\\b - X\\a) 2 
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For our value of the parameter a t , cos 0y ~ 1. Given the form (fT5|) . the intershell 
coupling (|16p becomes 

/•27T p2L a -(D{vi,v 2 )-A)/a t 

=W dvi dv 2 ^— 2 Ar e fai(gx>Jt+ax.fl.) e fa 3 (g Ht + gna ) (20a ) 

r 4n ,-2(L a +L b ) 

dui / d W2 e *M?^-^) e » 2 ( 9| ,a-^) =: A(q a , q b ) I(q a ,q b ), (206) 

Jo 

with vi = (xj_ a /R a - x±b/Rb), v 2 = x\\ a - x\\ b , and u% = (xj_ a /R a + x±_ b /R b ) and 
u 2 = i x \\a+ x \\b)- We denoted with A(q a , q b ) the amplitude of the coupling (it includes 
all numerical factors), while I(q a ,q b ) contains the functions which determine the 
selection rules discussed below and which appear upon performing the integration in 
(|206l) . It reads 

I(QaiQb) = 6 (K{q±aRa - q±bRb)) X 5 f La + Lb ( g||a _ g||6 )^ ; (21) 

where S(x) := sm(x)/x and q a = k a + G a ,q b = k b + G b . The resulting selection 
rules act differently on the angular and longitudinal degrees of freedom. It is when 
considering the latter that the issue of incommensurability arises. 
The additional integration over v\, v 2 yields the amplitude of the coupling A(q a , q b ). 
Although the integrals in (|20a[) are finite, the support of the integrand is well within 
the integration limits, which can therefore be extended to (— oo, oo). Thus we find 

A{q a ,q b ) = t k exp |~ ^^ (q± a R a + q±bRbf^ x exp |-^(g|| + <7|| b ) 2 j , (22) 

where tk ~ — 0.78eV contains both <o and all other numerical factors arising from 
the integrations. From (| 1 5 fl it follows that (|22[) and ([21]) have to be evaluated for 
q a = k a + G a , q b = k b + G b , with kp satisfying the boundary conditions (fT9]) . 
It clearly shows that contributions from distant regions of the momentum space are 
exponentially suppressed. For (k + G) > 27r/ao they are already negligible, therefore 
the sum in (|15p can be limited to only a few terms. 



3.1. Selection rules 

The selection function I(k a + G a , k b + G b ) determines whether the coupling between 
k a and k b is allowed. Note that all integer values of x/ir are zeroes of 5 = sm(x)/x, 
except x = where 6(0) — 1. 

Transverse degree of freedom. The angular momentum l_\_p — k±pRp can take only 
integer values and G±pRp = n\pm\p+n 2 pm 2 p G Z. Therefore S acts for the transverse 
degree of freedom q±R in the same way as a normal Dirac S. 

Longitudinal degree of freedom. The longitudinal momentum is k\\p = ^^l\\p- The 
longitudinal component of a reciprocal lattice vector is 

(Gp)n = (nxpbip + n 2 pb 2 p)\\ = Mp — \{2mxp + m 2 p)n 2 p - (m 1/3 + 2m 2 p)n 1 p)} , (23) 

L 

where Mp is the number of unit cells in shell (3, Lp — Mp\Tp\. The value of qup 
can therefore always be represented as ■j^d'p- If both shells are of equal lengths 
L a = L b = L, which is only possible in commensurate DWNTs, 

- ^La + Lb (na q ^ = ~ s{n(l ^ ^ = §{l ^ ^ (24) 
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In incommensurate DWNTs the two shells always have different lengths and the proper 
selection function is 5. However, as (|22l) shows, the amplitude of the coupling decreases 
strongly with the length of q a , q b , therefore only a finite region of the reciprocal space 
is active in the coupling, i.e. gives a non-vanishing contribution. If the mismatch 
between shell lengths is small enough, we can still approximate 5 by the Dirac S in 
the whole active region. Depending on which level of precision in this approximation 
we find acceptable, the maximum allowed mismatch can be large or small. We study 
only DWNTs for which (L a + Lb)(q\\ a — 9||b)/2 < 7r/4 in the active region and use the 
Dirac-5 selection rules below. 

l±a + (m la n la + m 2a n 2a ) = l±b + (mi b n lb + m 2b n 2b ), (25a) 

1 1| a + J 7 a (nia,n2a) = l\\b + Fb(nib, n 2 b) , (256) 

with 

Fpinip, n 2 p) = 2Mp ((2mip + m 2 p)nip - (my + 2m 2 p)n 2l3 ) . 

In the infinite DWNTs the integral over u 2 in (|206|) runs over infinity and 

the longitudinal selection rule becomes a true Dirac (5(q|| a — Q\\b)- The longitudinal 
momentum is continuous, which means that for any pair of momenta fc|| a , k\\ b exists at 
least one pair of G a , G b such that q^ a , q\\ b fulfill the selection rule, therefore in principle 
all longitudinal momentum states are coupled. However, in actual computation only 
the contributions from the active region of the reciprocal space count and many of the 
couplings vanish, reestablishing the division of the momentum space into independent 
subspaces, as was the case in finite DWNTs. 



3.2. Angular momentum: coupling between subbands 

In this section we analyze the implications of (|25a[) , i.e. we find the set of states 
which fulfill the angular momentum selection rules. As an example we consider the 
commensurate DWNT (5,5)@(10,10). The amplitude function A §FI§ is dominated by 
the regions in the reciprocal space corresponding to small values of G a ,G b . Among 
the reciprocal cells of dominant contributions wc find 

G a = (nio, n 2a ) = (0, 0), G b = (0, 0) l Xb = l± a 

ii) G„ = (-1,0) or (0,-1), G b = (-1,0) or (0,-1) =► lxb = l± a + 5 (26) 

ii) G a - (1, 0) or (0, 1), G b = (1, 0) or (0, 1) => l ±b = l ±a -5 

iv) G a = (1, 1) or (-1, -1), G b = (1, 1) or (-1,-1) => l± b = l± a - 10. 

For example, l± a = yields in case i) l±b = 0, in case ii) l± b = 5, in case Hi) l± b = —5 
and in case iv) l^ b = —10. These coupled states are shown in figure HI where for 
clarity only the states with h a = 0,lu b = are shown. Other combinations of G a ,G b 
would in this DWNT yield the same results. All combinations listed above fulfill also 
the second selection rule (|25b|) for l\\ a — l\\b- 

When we apply the selection rules (|25a|) in turn to all the found lj_ b , we find other 
l± a states which also couple to the lj_ b s found above. It turns out that in this 
particular DWNT the sets of coupled states contain only a few elements. For each 
initial l± a = l G [0,4], the set of coupled angular momentum values consists of 
l±a = — 5, Iq and l± b = Iq — 10, Iq — 5, Iq, Iq + 5. This is a rather unusual situation, 
occuring only when the chiral indices of one shell are integer multiples of those in the 
other. In an average DWNT the coupled sets are larger. 
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a: (5,0) -bj 

* a'.(5,0)-b 2 
b: (-5,0) + b 2 

° V: (-5,0) + b, 

c: (-10,0) + b, + b 2 
' c': (-10,0) 

• d: (0,0) 

Figure 4. Coupled subbands in a (5,5)@(10,10) DWNT. As an example we 
consider the point d = (l_i_a,l\\ a ) = (0, 0) a belonging to the reciprocal space of 
the inner shell. 61,62 are the graphene reciprocal lattice generators |(3j- After 
appropriate translations by graphene reciprocal lattice generators 61 , 62 ((3j the 
state (l± a ,l« a ) = (0, 0) a generates the inner shell states a,a',b,b' and c,c'. 
According to (1 26 I t each of these states couples to its counterpart a, ...,c' in the 
outer (10,10) shell, resulting in a non-zero coupling. 

3. 3. Longitudinal momentum - the issue of commensurability 

In our calculations we consider only DWNTs in which the two shells have equal or 
very similar length. If AL < oq/8 we are allowed to use <5-like selection rules and the 
physical space divides into subspaces containing the coupled longitudinal momentum 
states. The Hamiltonian matrix acquires a block-diagonal structure with the size of 
the blocks determined by the geometry of the shells. 

The active region of the reciprocal space can contain several reciprocal cells 
contributing to the coupling, which causes the mixing of longitudinal momentum 
states. The number of involved reciprocal cells and therefore of coupled longitudinal 
momenta increases with the size of the direct lattice unit cell. 

If the unit cell ratio of the two shells is rational, \ = vll where p, q £ N, the length of 
the DWNT is L = qM\T a \ = pM\Tf,\. The selection rules split the momentum space 
into M subspaces, each containing the full set of subbands ({l±}) for q longitudinal 
states in shell a and p longitudinal states in shell b. The size of each subspace is 
2{qS a +pS b ). 

In the case of incommensurate DWNTs, the lengths of the shells can be chosen so as 
to minimise AL and allow us to use the exact conservation of crystal momentum, 
as explained above. If these optimal values of shell lengths can be expressed as 
L a = qM\T a \ and Lf, — pM|Tf,|, where p, q € N, the Hamiltonian splits into M 
diagonal blocks, each of the size 2(qS a + pSt,)- The ratio p/q is in fact a rational 
approximation of the irrational x an d the precision of this approximation depends on 
the required value of the difference between shell lengths. 

The difference between the commensurate and incommensurate shells is shown in 
figure [5] in the case of two finite linear chains. 

3.4- Energy spectrum at the Fermi level 

The details of the spectrum at the Fermi level depend on the form of both intra- and 
intershell interaction, most notably on whether the curvature of the nanotube is taken 
into account or not. Among the effects of curvature in single- wall nanotubes are the 
rehybridization of a and tt bonds and varying angle between tt orbitals |44[ |45j . They 



(0,0) -b! 

■ (0,0) - b 2 
(0,0) + b 2 

• (0,0) + b t 
(0,0) + 

■ (0,0) -b r 
(0,0) 



(10,10) 
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Figure 5. Possible couplings between momentum states in commensurate and 
incommensurate finite chains. Solid black lines mark the first Brillouin zones, b a 
and bf, are the reciprocal lattice generators on chain a and b, respectively. Filled 
dots stand for states in the first Brillouin zones, open dots for the states in the 
rest of the reciprocal space. Dashed grey rectangles mark the regions in which the 
momenta on different chains match under <5. Black dashed arrows connect states 
which are coupled after the translational equivalence has been taken into account, 
a) Commensurate chains with unit cell ratio 2/3. The couplings between states 
l a , lb and 2 a ,2[, are the result of direct matching of momenta. The state 3 a has 
the same momentum as 3;,, but 3;, is equivalent to 1„ under the translation by b^, 
which means that 3 a , 1(, are also coupled. Similar situation occurs for 4 a and 4 B , 
which arc equivalent to l a and 2 ot and so forth. In the end, all states are coupled, 
although the coupling may be weak. It is possible to define a common Brillouin 
zone, with the length 2n/3a a = 2n/2a . b) Incommensurate chains with unit 
cell ratio l/\/3. The lengths of the chains cannot match - here they are chosen 
as 3 unit cells of chain a and 2 unit cells of chain b. The momenta in chain a are 
shifted with respect to those in chain b, but this mismatch is not large and the 
same couplings as in the case a) occur. 



result in variations in the bond length and bond angle between the lattice atoms, 
which can cause the opening of a diameter-dependent gap at the Fermi level in metal- 
lic SWNTs. Moreover, the band structure at the Fermi level depends on the relative 
position of the shells, as it was found to be the case in a (5,5)@(10,10) DWNT studied 
in [T7J [57]. When the (5,5)@10,10) DWNT is in a configuration of maximum sym- 
metry D§h [18l [46] . the only effect of the intershell hopping is a uniform split and 
shift of the Fermi subbands, resulting in the presence of four subband crossings. If 
the symmetry of the system is lowered, four pseudogaps (the largest of the order of 
0.1 eV) open in the spectrum [TT] . 

When, as in this work, curvature effects are neglected, only a subband shift is ob- 
served (see figure [8]) - in other words, our nanotube is always in the configuration of 
maximum symmetry. Due to the small size of the curvature-induced gap, we think 
that our model still yields a reliable description also of band features near the Fermi 
energy. 

The presence of a uniform shift between the subbands of a (5,5)@(10,10) DWNT can 
be understood by considering just the coupling between the Fermi subbands. The 
general Hamiltonian (fT4|) for k a and k;, becomes a 4x4 matrix if all other couplings 
are ignored. In the sublattice basis it has the form 



H(k a ,k b ) = 



( 
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where 7,3 (k^) = 7oX^=i exp(ik / 3 • dj) =: |— | exp(i0^), t a t is the coupling amplitude 
between k a ,k(, from (|22p and <^„„/ = iG a • r a „ — iGf, • t/,,,/ is the phase associated 
with hopping between different sublattices. It is clear that ipsB — Pab + ^Pba- When 
this Hamiltonian is expressed with the help of (jTTJ) in the valence/conduction basis, 
it becomes 

/ Yia\ f +H 

-|7-l ^-H 



T* 



f- 
f_ 



tab ( 

TV 

t a b 



matrix 



\fb\ 






-M / 



(28) 



The elements of the coupling matrix in this basis are 
t 



<PAb) I p i(VAB+VBA) 



~ 2 1 
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,i{e b -B a ) 



ab ( 

TV 



a i(9 b + ipBA) 



i{9b+VBA) 



-i(6 a —<pA. 



<) _|_ e K<PAB+<PBA)^ 
<) _|_ g«(¥>AB+¥>BA)^ 



(29) 
(30) 

(31) 
(32) 



In order to analyze the nature of the subband (anti)crossing at the Fermi level, we 
need to evaluate 7(k l a), t 6 5 <Pab and <^ba- Let us begin by the in-shell part. 
For Fermi subbands the angular momentum is set to l± a = —5 and l±b = —10. We can 
rewrite 7(k l a) as a function of the distance between longitudinal momentum and the 
Fermi point, Afc = fcii — kp. Both subbands have the same position in the reciprocal 
cell of the armchair nanotube, therefore 7 a (Afc) = 7k(Afc). The examination of 7(Afc) 
reveals that its phase has only two values: 



7 (Afc) = 



7o 



^2 cos f| + ^Afca J 



s i2 7r /3 ! Afc < 
s - 47 r/3 Afc > 



(33) 



In the inter-shell part we have to perform a sum over reciprocal lattice vectors as 
in (|15|) . The vectors which give the dominant contribution to the coupling are 
G a = bi,Gb = bi. The phases associated with hopping between sublattices are 
then 



*£AA = 0, IfiAB 

We have now a situation where 9 C 
matrix T becomes 



2tt 
T 

= b 



¥>BA 



2tt 

~T' 

6 and lpab = 



Vbb 
-<Pba 



0. 



(34) 

The coupling 



T = t r , 



1 + cos(6> — ip) 
i sm(8 — ip) 



-i sin(# — ip) 
— cos(8 — ip) 



= t n 







1(1 
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Us 
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(35) 



The coupling does not mix bands and moreover affects only the conduction band for 
k < kp and the valence band for fc > kp, as shown in figure [5] In consequence, the 
negative slope parts of the two subbands — 5 a and — IO5 are split evenly on both sides 
of the Fermi point. The results of our calculation show that there is also a smaller 
uniform split of the part with positive slope (see figure [8]), which is due to the smaller 



Electronic spectra of DWNTS in parallel magnetic field 



14 



conduction 



conduction 



a) 




b) 



T++ >0 



k 



= 



valence 



Figure 6. Schematic plot of the band crossing at the Fermi level in a 
(5,5)0(10,10) DWNT. a) The coupling between subbands at k < kf and k> k F . 
Only the conduction band is affected in the former, only the valence band in the 
latter case, b) Resulting shift of the energy levels close to the Fermi point. The 
parts with negative slope are split, the parts with positive slope remain degenerate. 



couplings between — 5 a and —5b and 5f>. This asymmetry is explained in [47j in real 
space terms as the result of different phases of the wavefunctions on sublattices A and 
B. The wavefunctions belonging to the negative slope parts of the subbands have 
constant phase on the whole circumference, while the wavefunctions belonging to the 
unshifted parts have different phases on sublattices A and B, therefore they cannot 
hybridize so well. 

3.5. Results 

In order to test our method we calculated the electronic spectra of a short commen- 
surate (5,5)0(10,10) and a short incommensurate (9,0)@(10,10) DWNT, both by the 
fc-space method described above and by direct diagonalization of the Hamiltonian in 
the real space. Our test DWNTs have both shells of equal or very similar length. The 
commensurate DWNT consists of 120 unit cell lengths of both the inner and outer 
armchair. The incommensurate DWNT has 75 unit cells of the zigzag shell and 130 
of the armchair. The mismatch between shell lengths is 0.17ao, which still allows us 
to use the exact selection rules ()25&[) . The spectra calculated by both methods match 
well (figure 0. 

The asymmetry between the valence (E < 0) and conduction (E > 0) bands, seen 
also in [25], is due to the intershell tunneling. The wave functions of the coupled 
momentum states hybridize and form bonding and antibonding combinations, with 
greatest energy differences in the bottom of the valence band. 

The coupling between momentum states is felt most strongly by states with low mo- 
mentum, but its consequences can also be seen at the Fermi level. In our model, which 
neglects curvature effects, the band structure at the Fermi level does not depend on 
the shell shift. In the k-space approach the relative position of the shells enters only 
through a phase factor in the coupling matrix (|15[) and does not affect the energy 
eigenvalues. We have tested our prediction for a (5,5)@(10,10) DWNT using a "par- 
tial real space" method. We defined a supercell containing one unit cell of the outer 
and one of the inner shell. Then we used the Bloch theorem in the longitudinal di- 
rection and the spectrum which we obtain is also insensitive to the shell shift (figure[5]). 
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Figure 7. Comparison between the DOS of a (5,5)@(10,10) (left panel) and 
a (9,0)®(10,10) (right panel) DWNT evaluated by diagonalizing the DWNT 
Hamiltonian in real (Eq. JEJl) and reciprocal (Eq. JHJl) space. n(E)/N 
is the density of states normalized to 1 (TV is the number of atoms in the 
nanotube). In (5,5)@(10,10) both shells have 120 unit cells and are 30 nm long. 
In the (9,0)@(10,10) the zigzag shell contains 75 unit cells, the armchair - 130, 
corresponding to the DWNT length of approximately 32 nm. Notice the breaking 
of the electron-hole symmetry due to the intershell tunneling. 




Figure 8. Anticrossing at the Fermi level in a (5,5)0(10,10) DWNT. For 
independent shells (a) the subbands at the Fermi level are degenerate, but when 
the intershell tunneling is allowed (b), the degeneracy is removed. The energy 
levels in the "partial real space" method were obtained by defining a DWNT 
supercell and using the Bloch theorem in the longitudinal direction. 



4. DWNT in parallel magnetic field 

The intershell coupling 

When a magnetic field is applied to a system, it usually changes the system's 
symmetries, since the vector potential A depends on the spatial coordinates. As 
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a consequence the wave function of a charged particle moving in the magnetic field 
gathers a phase factor during its motion. This is due to the modification of the 
momentum operator, p — > p — qA, known as the minimal coupling principle or the 
Peierls substitution. The translation operator T, where T(x)ip(r) = ip(r + x) (r is 
a position vector, x is the translation vector), is modified accordingly by the Peierls 
phase 28\: 

T(x)=exp^x.p) - r'(x)=cxp||^ r+X A(r')-dr'|r(x). (36) 

In a uniform field the Hamiltonian remains invariant under translations. For the 
lattices considered in the tight-binding model it implies that each hopping integral is 
modified by the appropriate phase factor. 

The influence of the magnetic field has been most extensively studied in two simplest 
cases - of a uniform magnetic field perpendicular to the flat lattice, and a uniform 
magnetic field parallel to the axis of a system with cylindrical topology. 
In the former case the application of the magnetic field changes or destroys the 
periodicity of a lattice. If the magnetic flux through the area of the elementary cell 
is rational, 4> ce ii = p/q <fto, where p, q € Z and (f>o = e/h is the flux quantum, it is 
possible to define an enlarged elementary cell, containing q original ones, pierced by 
p flux quanta. Thus the lattice remains periodic, although with a different period. 
If the flux through the elementary cell is irrational, 4>ceii/ ( t>o $ Qj the- periodicity is 
removed altogether and the spectrum of the system is fractal. The plots of the energy 
spectrum vs. 4> C eii are known as "Hofstadter butterflies" [29] . 

In systems such as rings and cylinders, a uniform magnetic field parallel to the axis 
gives rise to the Aharonov-Bohm effect. In simple systems its consequence is a shift 
of all the angular momentum states by the number of flux quanta flowing through its 
cross-section, 4> C ross- section I '4>o- The change in the spectrum is periodic with a period 
(f>o [33 132] • In DWNTs tunneling can also occur between shells and the Peierls 
phase enters not only into the in-shell term, but also into the inter-shell hopping (see 
figure [9|) and (fT2|) becomes 

H ( A ) = E E ^ exp j | I ~ A(r') • dr' \ c% a c^ 
P<y (id) 

+ 




(37) 

E *( r <»<' r bj) ex P \ V / A ( r ') ' dr ' \ c iia c bi<r + h.C. 



In the cylindrical coordinates and the tangential gauge we have chosen, A = 
(A r , A v , A z ) — (0, Br /2,0). The phase factors attached to in-shell bonds are 

{ie r T Pi 1 -*j3/ < j+P ( *±pi x ±nj \ 

j- J A(r')-dr'J = e l 7o(vpi-v») =e % 1 «i h„J ; ( 38 ) 

where <j>p is the magnetic flux through the shell (3. The dispersion relation contains 
therefore a dependence on the magnetic flux: 
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An 




Figure 9. DWNT in a uniform magnetic field parallel to its axis in tangential 
gauge. The atoms i a on the inner shell a and ji, on the outer shell b bind to 
their in-shell neighbours and with each other. A phase factor comes into the 
Hamiltonian with each bond. 



The phase factor in the term describing the inter-shell interaction is 
I " / " " . / , , I ' / -r±i,, .r±,„ \ I . Hi, i Hi, \ ' 

exp ■ 



h L, ^ ' J 3 00 V Rb Ra J \ Ra \Ra 



% _Lai 




(40) 



Rb 

The inter-shell coupling is analogous to (fT5|) except that t q q depends on the magnetic 
flux threading the DWNT. It is given by 

iq a ,q b (^a) = t k 5 (ir(q ±b R b - q± a Ra)) S ( La ~^ Lb (q\\ b - q\\ a )j (41a) 



Cxp { - ^ (i\\b + q\\a) 2 ^ 



(416) 



x ^\-^k{ q - bRb+ ^ A+2F {t)) }' (41c) 

where t% ~ 0.78eV. The magnetic field B enters the Hamiltonian only as the flux 
through inner and outer shell. The flux through the outer shell can be expressed as 
tfib = {Rb/ Ra) 2 4>a and in the following we shall present all quantities depending on the 
magnetic field as functions of the flux through the inner shell, <p a ■ 
As we can see from the 8 function in (141 ap . the selection rules do not depend on 
the magnetic field. They determine once and for all the quantum numbers of the 
coupled states, although for some pairs the main contribution to the sum in (|41a|l 
may come from a very distant reciprocal cell. The strength of the coupling, however, 
does depend on the amount of magnetic flux through the system. At vanishing field 
the most strongly coupled states are those with low momentum; at higher fields the 
maximum coupling can occur between states with energy close to the Fermi level or 
even to the top of the conduction band. From (|41cp we see that the strength of the 
coupling between angular momenta evolves with the magnetic field, while the coupling 
between longitudinal states (|416[) remains unchanged. It is to be expected, because 
the longitudinal motion of the electron does not accumulate the Peierls phase. It is 
clear that as the flux through the DWNT is increased, the dominant terms in the sum 
(|41ap come from reciprocal cells with varying G± but constant G\\ = 0. 
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Figure 10. Subbands coupled to lj_ a = in the inner shell of a (5,5)0(10,10) 
DWNT and the reciprocal cells of the main contributions. As the magnetic field 
increases, more distant cells are involved and the dominant coupling switches 
between (l± a ,i±b) = (0a,0(,) and (0 a , — 10b) with the period <E> w 4.29</>o- 

Let us analyse the influence of the magnetic field on the coupling between individual 
subbands. As an example we take the l± a — and l± b = — 10, Zj_{> = subbands of 
the (5,5)@(10,10) DWNT. The value of l\\ a = l\\ b = 1 shall be assumed implicitly. In 
the absence of the magnetic field the coupling between subbands (0 a , Of,) is dominant, 
while the coupling between (0 a ,— 10b) almost vanishes. As we increase the magnetic 
field, the predominant coupling switches between (0 a ,0b) and (0 Q ,— 10b), while also 
oscillating in amplitude. The switching occurs periodically and the period can be 
evaluated from (|41dl . The maxima of the coupling occur when the exponent vanishes: 

(l± a + n a S a ) + (l ±b + n b S b ) + + ^ + ^ = 0. (42) 

The first maximum of (0 o , b ) coupling occurs at (f> a = 0. The next maximum coupling 
is between (0 o , — 10^) and occurs at 

-10-10 + 4.67%^ ^</Wx« 4.2900- 

The period of the oscillation of the coupling amplitude is <E> « 4.290o- The switching 
between dominant couplings ((0 a ,0t,) and (0 a ,— 10&)), depending on which reciprocal 
cell is active, is shown in figure [TO) With the magnetic field increasing from 0, the 
reciprocal cells of the dominant contribution change in a sequence 

<t>a =0, G a = 0, G b = 

a = $, G a = -l(6i + b 2 ), G b = 

<p a = 2n$, G a = -2n(& a + b 2 ), G b = + b 2 ) 

4> a = (2n + 1)$, G a = -(2n + l)(6i + b 2 ), G b = -n(bi + b 2 ), 

where n — 1,2.... If the origins of the shells are aligned, p a = p b = (see figure [3]), 
the phase factors from (TT5)) change with the period 6$, common to all pairs of coupled 
states. If any of the p's is non-zero, the factor exp(zG • p) is periodic in cf> only for 
p = qao, with q rational. Otherwise the phase factors in (|15p vary in the magnetic 
field without showing any periodicity. 

In DWNTs where the chiral indices of the outer shell are not integer multiples of those 
in the inner shell, the coupling between lj_ a = and its partners in the outer shell 
also strengthens and weakens periodically. For example, in a (6,6)@(11,11) armchair 
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nanotube the subband a is coupled to all the even-numbered subbands in the outer 
shell. In the absence of the magnetic field the dominant coupled pair is (0 a , Ob). As the 
field increases, the dominant pair becomes (0 a , 10b), then (0 a , — 2&), and so forth. The 
distance between subsequent maxima of the coupling strength can be found as above, 
by minimizing the exponent in (|41 c|) and in this particular DWNT it is $ w 5.810o- In 
an incommensurate (9,0)@(10,10) the situation is analogous and the oscillation period 
is 8.1600- 

4-2. Results 

We performed numerical calculations of the DWNT spectra and explored the evolution 
of coupled states for several combinations of chiralities. In figure [Tl] we show the 
behaviour of the coupling between two sets of states of a (5,5)@(10,10) in the 
magnetic field. Those sets are (l a ±,l a \\) G {(— 5, 1), (0, 1)} in the inner shell and 
{h±,lb\\) £ {(— 10, 1), (— 5, 1), (0, 1), (5, 1)} in the outer. In the regions where the 
coupled states have both similar energies and strong coupling we notice the appearance 
of avoided crossings. Their size and position in the spectrum is governed by four factors 
depending on the magnetic field: 

i) the amplitude of the coupling (|415[) . (141c[) 

ii) the phase factors (|T5|) 

Hi) the dispersion relations e a ^(k, a /</>o) 
iv) and eb,i>{k,<j> a / 4>o) (Ell- 
in this particular case the spectrum is periodic in cj) a - the strength of the coupling 
oscillates with the period $ = 30(/>o/7, the phase factors with the period 6$, the energy 
of the inner states with the period <pOi the energy of the outer states with the period 
(Rb/R a ) 2 4>o — 40o- At cj) a — 18000 the initial spectrum is recovered. Nevertheless, 
in other DWNTs Rb / R a is usually irrational and it is in general impossible to find a 
common period for these four functions. 

The numerically calculated DOS plots of several nanotubes show features absent in 
uncoupled DWNTs (see figure [P2|) . The coupling between states from the two shells 
causes a series of avoided crossings, resulting in a whole region in which the density 
of states is depleted, observed also in [3S]. In small magnetic field this region is at 
the bottom of the valence band, where the momenta in both shells are small. As the 
inter-shell tunneling evolves with the increasing magnetic field and higher momentum 
states become involved, the main region with avoided crossings shifts also towards 
higher energies. 

These snake-like patterns are a statistical result, caused by many states. Their details 
vary according to the chiralities of the DWNT's shells. In armchair nanotubes the 
most strongly coupled states are at the band edges, and the avoided crossing affects 
van Hove singularities. In other nanotubes the strongest coupling can occur farther 
from the band edges, especially when the magnetic field is large. These energy gaps 
evolving in the middle of the band give rise to less distinctive snake-structures which 
can be seen seen in figure 1131 where they are caused by the presence of a zigzag shell. 

The evolution of the density of states with the magnetic field is a superposition 
of two patterns at different scales. Features with steep E/<fi slope are caused by the 
outer shell, which feels a flux greater than that in the inner shell and evolves faster 
with the magnetic field. The features with mild slope are due to the inner shell. 
Characteristic of the evolution of the DOS near the Fermi level with the magnetic field 
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Figure 11. The evolution of two sets of coupled states in a (5,5)@(10,10) 
DWNT. The inner shell states are (lxaJ\\a) 6 {(-5, 1), (0, 1)}, the outer shell 
Qxb> hlb) E {( — 10, 1), (—5, 1), (0, 1), (5, 1)}. The abscissa corresponds to the inner 
flux in <po units, the ordinate to the energy, a) The energy of the states (—5, 1);, 
(black), (5, 1);, (light grey) and (—5, l) a - The latter is shown either in dark grey, if 
the prevailing coupling is (—5, l) a with (—5, l)t or in medium grey if (—5, l) a with 
(5, l);, dominates. The width of the lines gives additional information about the 
size of the coupling, e.g. at </> a « 240o the prevailing coupling is between (—5, l) a 
and (5, 1);, (wide medium grey and light grey lines), which are also close in energy. 
The state (—5, 1){, has a much higher energy and couples to (—5, l) a very weakly 
(small black dots). At border values of <j> a where the dominant coupling switches 
phase, black lines are drawn, b) The difference between the energy spectrum 
for this subspace obtained without inter-shell tunneling (to = 0) and with the 
tunneling of the magnitude to = 70/8. c) The analogon of a), but for the sets 
of states (—10, 1);, (black), (0, 1)(, (light grey) and (0, l) a (medium or dark grey, 
depending on the prevalent coupling) . It can be seen that at flux values where the 
energies of strongly coupled states are close, there is a distinctive avoided crossing 
in the energy spectrum, such as e.g. at <j> a = or <j> a ?s 4</>n in c) or 4> a S3 24</>o 
and <f> a ~ 360o in a). The corresponding regions of large avoided crossings arc 
marked. Whether the crossing occurs in valence or conduction band is determined 
by the phase factor II16I I for the Brillouin zone of the dominant contribution. 
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Figure 12. The density of states in a (6,6)0(10,10) DWNT in changing parallel 
magnetic field, a) Without inter-shell tunneling the DOS is a sum of the DOS in 
both shells, b) The inter-shell coupling causes a change in the DOS varying with 
the strength of the magnetic field. Here the coupling constant is to = 7/8. 
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Figure 13. The DOS of an incommensurate (9,0)0(10,10) DWNT in changing 
magnetic field. Note that in most of the energy range the density of states seems 
decreased in comparison with the commensurate case in Fig. 1121 That effect is 
due to a very high DOS at the van Hove singularities occurring at ±7, typical for 
the zigzag nanotubes. 

is the periodic opening and closing of the gap, causing a series of metal-semiconductor 
transitions, predicted in [48] and observed a few years ago in [33l 134] . These can be 
seen as the empty diamonds along the E = Q line in figure [T4l In a (6,6)@(11,11) 
nanotube the intershell tunneling mixes the subbands from both shells, increasing the 
DOS at the Fermi level whenever both nanotubes have closed gaps, e.g. at <j) = and 
4> w 20o in figure [T4"a). In an incommensurate (10,0)@(11,11) DWNT the coupling 
between shells affects at the Fermi level only the DOS of the inner shell - due to the 
coupling the band crossing is shifted towards positive energies - the large diamonds 
are shifted with respect to the small ones in figure [14b) . This is a consequence of the 
structure of the reciprocal cell of the zigzag shell, where the Fermi point (subbands 
crossing) is at In = 0, while for the armchair it is at l\\ = 21f/3. The amplitude of 
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the coupling decreases exponentially with the value of the momentum, therefore the 
zigzag subbands at the Fermi level are affected more strongly by the coupling than the 
armchair subbands. In transport experiments, where mostly the outer shell is probed, 
this shift might be visible if the zigzag shell is on the outside. 



(6,6)18(11,11) ' (10,0)@(11, 11) 




1 * ,a 2 1 2 

<M4>o 4>„/0 o 



Figure 14. The evolution of the DOS at the Fermi level with magnetic field in a 
commensurate (6,6)0(11,11) (a) and incommensurate (10,0)0(11,11) (b) DWNT. 
The closing and opening of the gap can be seen in both cases. In (a) the peaks 
belonging to different shells are mixed, in (b) the only effect of the intershell 
coupling is a shift of the band crossing in the zigzag shell. 



5. Conclusions 

In this work we started from an inter-shell tunneling Hamiltonian given in the 
real space and derived its equivalent in the reciprocal space. In a commensurate 
(5,5)@(10,10) DWNT the band structure obtained with this method agrees with that 
obtained by the partial real-space method described in section 13.51 down to the fine 
details of the subband crossings near the Fermi level. As shown in section I3.4| this 
method allows us also to study the spectrum near the Fermi level analytically. Al- 
though for small nanotubes the curvature (which we neglect) can cause a dependence 
of the spectrum at the Fermi level on the relative position of the shells [17] , we expect 
this effect to decrease strongly with the nanotube diameter. Our method is therefore 
suitable for the realistic DWNTs with diameters above 2nm [39] . 
When this method is applied to the DWNTs in the parallel magnetic field, we observe 
complex geometrical patterns developing in the DOS of the nanotubes. The most 
prominent ones are at energies inaccessible experimentally, but we find the effects of 
the intershell coupling also at the Fermi level. In a double-armchair DWNT we find 
the metallic character of the tube persisting also at <f> > 0, while without this coupling 
the system would become semiconducting immediately after switching on the field. In 
an incommensurate zigzag@armchair DWNT we find the band crossing of the zigzag 
shell shifted towards higher energies, while the band structure at the Fermi level in 
the armchair shell is almost unaffected by the intershell coupling. 
The real-space methods of finding the spectrum of long commensurate DWNTs, where 
it is possible to define a common unit cell, are usually fast enough. In the case of in- 
commensurate DWNTs the real space approach must be either to diagonalize the 
Hamiltonian of the whole DWNT, or to squeeze or stretch one of the shells so that 
they become commensurate and an approximate supercell can be found. The former 
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is very costly in terms of computation time and memory, the second involves a de- 
formation of the lattice of one or both shells. Solving the Schrodinger equation in 
momentum space, as described here, allows us to use the selection rules and signifi- 
cantly reduce the size of the matrices to diagonalize. This method has been proven 
correct for short incommensurate nanotubes and for long commensurate DWNTs in 
a parallel magnetic field, where it gives the same results as those obtained in [39]. It 
may be a useful tool in investigating other properties of DWNTs. 
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Appendix A. Derivation of the intershell coupling matrix elements 

In this appendix we give an explicit derivation of the elements of the intershell coupling 
matrix, (|15[) . By expressing the Hamiltonian in the basis of plane waves in each shell, 
the intershell coupling Hamiltonian H t between two shells a and b can be written as 

Ht =^U ai ,r bj c\ ia c bj(T + H.c. 

ijcr 

— 1 . \^ \ " -ik-r ai +ik-r bj , t , it 

~~ . 'V A T ' ' ' < L r az ,r b:I ^ apaka< jHp b k b o T xi.c, 

V^a^b ijcr k a k bPaPb 

where p a ,Pb = ± are the indices for the two interpenetrating sublattices, N a and N b 
are the number of graphene unit cells on shell a and 6, respectively. A carbon atom 
can be found in graphene at the position r = R+ p+pr, where R is a graphene lattice 
vector and p and r are two vectors that specify the atom position, cf . figure [3Ja) . 
The sum over all the lattice sites can be carried out as 

^ \ " ~ik-r ai +ik-r b jj. 



^^R a R bPaPb 



^ ' e -ik-r ai +ik-r bj ^ 



r a =R a +p+ Pa T,r b =R b +p+ Pb T 



==X] / dr a dr b $( r a- (Ra+ P+PaT)) 

a-^b _ „. J d d 



x S(r b -(R b + p + p b T))e- lk - r " +lk ^t rat , rbj 

_ 1 e iG a -(p+p a T-)- l G 6 -(p+p () -r) 



r ai ,r bi i 



where we have replaced the sum over the lattice vectors R by the sum over the 
reciprocal lattice vector G as 

£ 5(r (R + p + P t)) = J- £ e-iCCr-p-pr) 

R ^ coU G 

with .i4. ce n the area of a graphene unit cell. Thus we obtain the expression of the 
elements of the intershell coupling matrix, f) 15[) . 



